The Infrared Behavior of Gluon and Ghost Propagators in Landau Gauge QCD 
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A solvable systematic truncation scheme for the Dyson-Schwinger equations of Euclidean QCD 
in Landau gauge is presented. It implements the Slavnov-Taylor identities for the three-gluon 
and ghost-gluon vertices, whereas irreducible four-gluon couplings as well as the gluon-ghost and 
ghost-ghost scattering kernels are neglected. The infrared behavior of gluon and ghost propagators 
is obtained analytically: The gluon propagator vanishes for small spacelike momenta whereas the 
ghost propagator diverges more strongly than a massless particle pole. The numerical solutions are 
compared with recent lattice data for these propagators. The running coupling of the renormalization 
scheme approaches a fixed point, a c ~ 9.5, in the infrared. 
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A theoretical understanding of confinement of quarks 
and gluons into colorless hadrons could be obtained by 
proving the failure of the cluster decomposition property 
for color-nonsinglet gauge-covariant operators. One long 
established idea in this direction is based on the occur- 
rence of infrared divergences to suppress the emission of 
colored states from color-singlet states [Q. Such a de- 
scription of confinement in terms of perturbation theory 
necessarily has to fail. 

Thus, to study the infrared behavior of QCD ampli- 
tudes non-perturbative methods are required, and, since 
divergences are anticipated, a formulation in the contin- 
uum is desirable. Both of these are provided by stud- 
ies of truncated systems of Dyson-Schwinger equations 
(DSEs), the equations of motion of QCD Green's func- 
tions. Typically, for their truncation, additional sources 
of information like the Slavnov-Taylor identities, entailed 
by gauge invariance, are used to express vertex functions 
in terms of the elementary 2-point functions, i.e., the 
quark, ghost and gluon propagators. Those propaga- 
tors can then be obtained as selfconsistent solutions to 
non-linear integral equations representing a closed set 
of truncated DSEs. Some systematic control over the 
truncating assumptions can be obtained by successively 
including higher n-point functions in selfconsistent calcu- 
lations, and by assessing their influence on lower n-point 
functions in this way. At present, even at the level of 
propagators no complete solution to truncated DSEs of 
QCD exists. In particular, even in absence of quarks, so- 
lutions for the gluon propagator in Landau gauge rely on 
neglecting ghost contributions J^-^]. While this particu- 
lar problem is avoided in ghost free gauges such as the 
axial gauge, in studies of the gluon DSE in this gauge ||, 
the possible occurrence of an independent second term in 
the tensor structure of the gluon propagator has so far 
been disregarded jf) . In fact, if the complete tensor struc- 
ture of the gluon propagator in axial gauge is taken into 
account, one arrives at equations of no less complexity 
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than the ghost-gluon system in the Landau gauge || . 

In addition to the prospect of some insight into con- 
finement from studying the infrared behavior of QCD 
Green's functions, DSEs have proved to be a highly 
successful tool in developing a hadron phenomenology 
that interpolates smoothly between the infrared (non- 
perturbative) and ultraviolet (perturbative) regimes ||. 
In particular, a variety of models for the interactions of 
quarks mediated by gluons exists, which are very well 
suited for a dynamical description of chiral symmetry 
breaking from the DSE of the quark propagator in some 
analogy to the gap equation in superconductivity [ fl0| |. 
The superficial result of these studies is that for the 
quark self-energy to reflect a spontaneous breaking of 
chiral symmetry there has to be some sufficient interac- 
tion strength at low energies. Under these circumstances, 
the dichotomy of the pion as a Goldstone boson emerging 
from the Bethe-Salpeter equation for quark-antiquark 
bound states is very well understood and explains the 
smallness of its mass as compared to all other hadrons. 

In this letter we present a simultaneous solution of 
a truncated set of DSEs for the propagators of gluons 
and ghosts in Landau gauge. An extension to this self- 
consistent framework to include quarks dynamically is 
possible and subject to further studies. The behavior 
of the solutions in the infrared, implying the existence 
of a fixed point at a critical coupling a c ~ 9.5, is ob- 
tained analytically. The gluon propagator is shown to 
vanish for small spacelike momenta in the present trun- 
cation scheme. This behavior, though in contradiction 
with many previous DSE studies 11, can be partially 
understood from the observation that, in our present cal- 
culation, the previously neglected ghost propagator as- 
sumes an infrared enhancement similar to what was then 
obtained for the gluon. 

Besides all elementary 2-point functions, i.e., the 
quark, ghost and gluon propagators, the DSE for the 
gluon propagator also involves the 3- and 4-point ver- 



tex functions which obey their own DSEs. These equa- 
tions involve successively higher n-point functions. A 
first step towards a truncation of the gluon equation is 
to neglect all terms with 4-gluon vertices. These are 
the momentum independent tadpole term, an irrelevant 
constant which vanishes perturbatively in Landau gauge, 
and explicit 2-loop contributions to the gluon DSE. The 
latter are subdominant in the ultraviolet and will thus 
not affect the behavior of the solutions for asymptoti- 
cally high momenta. In the infrared it has been argued 
that the singularity structure of the 2-loop terms does 
not interfere with the one-loop terms [Q . Without con- 
tributions from 4-gluon vertices (and quarks) the renor- 
malized equation for the inverse gluon propagator in Eu- 
clidean momentum space is given by |l3| ], 
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where p = k + q, D tl and T tl are the tree level propagator 
and 3-gluon vertex, D G is the ghost propagator and T 
and G are the fully dressed 3-point vertex functions. The 
equation for the ghost propagator in Landau gauge QCD, 
without any truncations, is given by 
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The renormalized propagators for ghosts and gluons and 
the renormalized coupling are defined from the respec- 
tive bare quantities byjntroducing multiplicative renor- 

Z 3 D^ := D% and 
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The ghost 

and gluon propagators are parameterized by their respec- 
tive renormalization functions G and Z , 



malization constants, Z%D G 
Z g g := go. Furthermore, Z\ = Z g Z^'^, Z\ 
and we use that Z\ = \'va Landau gauge |l 
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In order to arrive at a closed set of equations for the func- 



tions G and Z, we use a form for the ghost-gluon vertex 
which is based on a construction from its Slavnov-Taylor 
identity (STI) neglecting irreducible 4-ghost correlations 
in agreement with the present level of truncation [fl5f , 
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With this result, we can construct the 3-gluon vertex ac- 
cording to general procedures from previous studies p6[ , 



wp{p,q,k) = ^A+(p 2 ,q 2 



k 2 ) V i(p - q) f 



(5) 



1 



A_(p 2 ,q 2 ;k 2 



2 

x (ti^pq 



)S^i{p + 

Pvq^)i{p- q) P 



q) f 



A_(p 2 ,q 2 ;k 2 ) 



cyclic permutations , 



A ± {p 2 ,q 2 ;k 2 ) 



G{k 2 )G{g 2 ) G(k 2 )G(p 2 ) 
G(p 2 )Z(p 2 ) G(q 2 )Z(q 2 ) 



Some additionally possible terms, transverse with respect 
to all three gluon momenta, cannot be constrained by its 
STI and are thus disregarded. For the fermion vertex in 
QED as constructed from its Ward-Takahashi identity 
it is well known that additional transverse terms, with 
the further constraint not to introduce kinematic singu- 
larities, are essential for multiplicative renormalizability 
]l7t . Based on this requirement such terms have been ob- 
tained explicitly for quenched QED in ref. p8| . Similar 
constructions for the vertices in QCD are presently not 
available. However, the full Bose (exchange) symmetry 
of the 3-gluon vertex alleviates this problem since, com- 
bined with the STI, it puts much tighter constraints on 
this vertex then those obtained for fermion vertices. 

Instead of a direct numerical solution of the coupled 
system of integral equations resulting from the present 
truncation scheme we use a one-dimensional approxima- 
tion: For integration momenta q 2 < k 2 we use the an- 
gle approximation replacing G((k — q) 2 ) — ► G(k 2 ) and 
Z((k — q) 2 ) — ► Z(k 2 ). Since this preserves the limit 
q 2 — > 0, it is suitable for an analytic discussion of the 
solutions in the infrared. For q 2 > k 2 we replace all ar- 
guments (including the external k 2 ) by the integration 
momentum q 2 . The justification for this is the weak log- 
arithmic momentum dependence of G and Z at high mo- 
menta @. The DSEs (§) and (|) then simplify to 
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We introduced an 0(4)-invariant momentum cutoff 
A uv to account for logarithmic ultraviolet divergences 
which are absorbed by the renormalization constants Z% 
and Z$. Z\ has to be ultraviolet finite 0). This is in- 
consistent with gauge invariance implying Z\ = Z3/Z3. 
While this problem, appearing at order g in a perturba- 
tive expansion, is quite natural for a truncation scheme 
neglecting explicit 4-gluon couplings at the same order, 
its remedy could provide information on purely trans- 
verse terms in the 3-gluon vertex. For details of the 
renormalization and the numerical procedure see p(jf . 

To deduce the infrared behavior of the propagators 
we make the Ansatz that for x := k 2 — > the prod- 
uct Z(x)G(x) — * cx K with k 7^ and some constant c. 
The special case k — leads to a logarithmic singularity 
in eq. (fjj) for x — ► which precludes the possibility of 
a selfconsistent solution. In order to obtain a positive 
definite function G(x) for positive x from an equally pos- 
itive Z{x), as x — > 0, we obtain the further restriction 
< k < 2 . Eq. (0) then yields, 

G(x)^( 9 * 7 $(--l)Y 1 c-ix-« =► (8) 
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where 7^ = 9/(647r 2 ) is the leading perturbative coef- 
ficient of the anomalous dimension of the ghost field. 
Using (W) and (0) in eq. (||), we find that the 3-gluon 
loop contributes terms ~ x K to the gluon equation for 
x — > while the dominant (infrared singular) contribu- 
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Comparing this to (^|) we obtain a quadratic equation 
with a unique solution k — (61 — Vl897)/19 ~ 0.92 for 
the exponent k < 2 . The leading behavior of the gluon 
and ghost renormalization functions is entirely due to 
ghost contributions. The details of the approximations 
to the 3-gluon loop have no influence on these consid- 
erations. In particular, additional transverse terms of 
the 3-gluon vertex, free of kinematical singularities, will 
yield contributions that are even further suppressed in 
the infrared. Compared to the Mandelstam approxima- 
tion, in which the 3-gluon loop alone determines the in- 
frared behavior of the gluon propagator and the running 
coupling in Landau gauge [^|-^), this shows the impor- 
tance of ghosts. The result presented here implies an 
infrared stable fixed point in the non-perturbative run- 
ning coupling of our subtraction scheme, defined by 



a s (s) = -Z(s)G^s) ^ — --- 



9.5 , (10) 



for s — » 0. This is qualitatively different from the infrared 
singular coupling of the Mandelstam approximation H . 



The momentum scale in our calculations is fixed from 
the phenomenological value as(Mz) — 0.118 at the mass 
of the Z-boson |2l}| . The ratio of the Z— to the r-mass, 
M Z /M T ~ 51.5, then yields a s (M T ) = 0.38 which is in 
encouraging agreement with the experimental value. 

It is interesting to compare our solutions to recent lat- 
tice results using implementations of the Landau gauge 
condition §| || . In figure Q we compare our solution for 
the gluon propagator to data from ref. p3] . We normal- 
ized the gluon propagator according to Z(x = 1) » 11.3 
to account for the units used in ref. (with x — k 2 a 2 
in units of the inverse lattice spacing) . According to the 
authors of [23[ , the arrow indicates a bound below which 
finite size effects become considerable. 

In figure |^ we compared our infrared enhanced ghost 
propagator to the results of ref. p3] . It is quite amazing 
to observe that our solution fits the lattice data at low 
momenta significantly better than the fit to an infrared 
singular form Dc{k 2 ) = c/k 2 + d/k 4 given in Q. We 
therefore conclude that present lattice calculations con- 
firm the existence of an infrared enhanced ghost prop- 
agator of the form D G ~ l/(k 2 ) 1+K with < k < 1. 
This is an interesting result for yet another reason: In 
ref. |Q the Landau gauge condition was supplemented 
by an algorithm to select gauge field configurations from 
the fundamental modular region which is to eliminate 
systematic errors that might occur due to the presence 
of Gribov copies. Thus, our results suggest that the ex- 
istence of such copies of gauge configurations might have 
little effect on the solutions to Landau gauge DSEs |25|| . 

The Euclidean gluon correlation function presented 
here can be shown to violate reflection positivity p0[ , 
which is a necessary and sufficient condition for the exis- 
tence of a Lehmann representation |^6| . We interpret this 
as representing confined gluons. In order to understand 
how these correlations can give rise to confinement of 
quarks, it will be necessary to include the quark propaga- 
tor. The size of the coupling at the fixed point, a c 9.5, 
is however, a good indication that dynamical chiral sym- 
metry breaking will be generated in the quark DSE. 

In summary, we presented a solution to a truncated 
set of coupled Dyson-Schwinger equations for gluons and 
ghosts in Landau gauge. The infrared behavior of this 
solution, obtained analytically, represents a strongly in- 
frared enhanced ghost propagator and an infrared vanish- 
ing gluon propagator. Our results, in particular for the 
ghost propagator, compare favorably with recent lattice 
calculations |^3|,^| . Since the lattice implementations of 
the Landau gauge are such that configurations are re- 
stricted to the fundamental modular region, this might 
indicate that Gribov copies have little influence on solu- 
tions to the DSEs in Landau gauge. The absence of a 
Lehmann representation for the gluon propagator can be 
interpreted as signal for confined gluons. The existence 
of an infrared fixed point is in qualitative disagreement 
with previous studies of the gluon DSE neglecting ghost 



contributions in Landau gauge This shows that 

ghosts are important, in particular, at low energy scales 
relevant to hadronic observables. 
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FIG. 1. The numerical result for the gluon propagator from 
Dyson-Schwinger equations (solid line) compared to lattice 
data from fig. 3 in |23| . 
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FIG. 2. The numerical result for the ghost propagator from 
Dyson-Schwinger equations (solid line) compared to data 
from fig. 1 in [^| for the 24 4 lattice up to i ~ 1, and a 
fit as obtained in ref. |24| for x > 2. 



